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The superconductor-to-insulator quantum phase transition in resistively shunted Josephson junc- 
tions is investigated by means of path-integral Monte Carlo simulations. Our numerical technique 
allows to directly access the (previously unexplored) regime of the Josephson-to-charging energy 
ratios Ej/Ec of order one. Our results unambiguously support an earlier theoretical conjecture, 
based on renormalization-group calculations, that at T — > the dissipative phase transition occurs 
at a universal value of the shunt resistance Rs = h/4e 2 for all values Ej/Ec- On the other hand, 
finite-temperature effects are shown to turn this phase transition into a crossover, which position 
depends significantly on Ej/Ec, as well as on the dissipation strength and on temperature. The 
latter effect needs to be taken into account in order to reconcile earlier theoretical predictions with 
recent experimental results. 



I. INTRODUCTION 

Mesoscopic Josephson junctions arc well known to ex- 
hibit a variety of intriguing phenomena which are of pri- 
mary importance both from a fundamental point of view, 
as well as for various applications including quantum- 
state engineering with electronic devices. Among these 
are macroscopic quantum tunneling with dissipation,!!] 
Coulomb blockade pa macroscopic quantum coherence, 
and dissipative quantum phase transitions. a Recent 
progress in nanolithographic techniques allows one to 
routinely fabricate ultrasmall tunnel junctions with ca- 
pacitances C as low as 10 -15 — 10~ 16 F, and to perform 
detailed experimental studies of various features related 
to the above phenomena. 

In the course of these studies it was realized - both the- 
oretically and experimentally - that the observed prop- 
erties of the system may crucially depend on the na- 
ture of the effective electromagnetic environment coupled 
to a mesoscopic junction. One of the most remarkable 
consequences of this dependence for Josephson junctions 
(JJ) is the quantum (T — 0) superconductor-to-insulator 
phase transition (SIT) driven by dissipation. The latter 
is controlled, e.g., by the magnitude of the ohmic shunt 
resistance Rs of the external leads. "Shis quantum phase 
transition was predicted by SchmioH and subsequently 
studied in Rcfs. (see also Refs. ||] for an extensive 
review of this and later theoretical activity). Thus, at low 
temperatures the supercurrent in mesoscopic JJs can be 
maintained only provided that quantum fluctuations of 
the Josephson phase are suppressed by dissipation. If 
dissipation is not strong enough, quantum fluctuations 
wash out the Josephson effect and no supercurrent can 
flow across the system. Alternatively, one can interpret 
this SIT as a destruction of Coulomb blockade for Cooper 



pairs by quantum fluctuations of the charge in an exter- 
nal resistor Rs- Experimentally this quantum dissipa- 
tive phase transition for single resistively shunted JJs was 
studied in Refs. 10- 13. The results of these experiments 
are qualitatively consistent with the above physical pic- 
ture. 

Similar experimental studies have been also performed 
for J J arrays and chains oli3 In that case an interplay 
between short and long-range quantum fluctuations of 
the superconducting phase in the presence of dissipation, 
yields a nontrivial low-temperature phase diagram.EZrt£l 
A quantum dissipative phase transition was also dis- 
cussed in the . casci of ultra-thin homogeneous supercon- 
ducting wires 

It is worth to point out that the above physical picture 
is not restricted to superconducting systems only. For in- 
stance, it is well known that the problem of a quantum 
resistively shunted JJ is equivalent to that of a quan- 
tum particle diffusing in a periodic potential coupled to 
a dissipative environment. In this case, the phase transi- 
tion from diffusion to localization occurs upon increasing 
the carupling strength to a dissipative Caldeira-Leggett 
bath.ErEI A formally identical Lagrangian also describes 
tunneling of electrons in a Luttinger liquid; see, e.g., 
Ref. 0. Similar physics was discussed for normal metallic 
conductors .E3 Thus, even though below we will specif- 
ically address the case of a resistively shunted J J, our 
results can also be applied in other physical situations. 

According to the existing theoretical picturejli at T = 
quantum localization of the Josephson phase should oc- 
cur as the shunt resistance becomes equal to the quantum 
resistance unit, Rs — R q = irh/2e 2 ~ 6.5 kf2, indepen- 
dently of the strength of the Josephson coupling Ej. In 
the limit of small Josephson energies (as compared to the 
effective charging energy of the junction, Ec — e 2 /2C) 
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this conclusion can be justified within a perturbative 
renormalization group analysis. Such an analysis can 
then be extended to the limit of large Ej ^> Ec by 
means of a duality transformation between the phase 
and the chargeBcl Since in both limits one obtains an 
-Ej-independent phase boundary at Rs = R q , it is rea- 
sonable to conjecture that its position remains unchanged 
for all values Ej/Eq including the regime of experimen- 
tal interest Ej ~ Ec- 

Although this conjecture can further be supported by 
a number of qualitative arguments, as well as by the 
existence of a self-duality point the structure of the 
phase diagram for moderate values Ej/Ec still needs 
to be rigorously, verified. Moreover, the results of re- 



cent experimentsEflliij could be interpreted as contradict- 
ing the above conjecture. In these works, the samples 
with sufficiently large Ej /Ec ^ 7 -f - 8 were found to be 
superconducting even for shunt resistances Rs substan- 
tially higher than R q . This observation could suggest 
that the true phase boundary should depend not only 
on the amount of dissipation in the system, but also on 
the ratio Ej/Ec- A similar conclusion could be reached 
from the results reported in Refs. |Io| , |l5| . 

All these developments motivated us to perform an ad- 
ditional theoretical investigation of the dissipative phase 
transition in a single resistively-shunted superconducting 
junction, at moderate values of the Josephson coupling 
energy Ej ~ Ec- Since in this range there exists no small 
parameter in the problem, it can hardly be rigorously 
investigated by analytical methods. Therefore, in this 
paper we analyze the problem numerically by means of 
path-integral (PI) Monte Carlo (MC) simulations. This 
method allows us to quantitatively study the effect of 
quantum fluctuations of the Josephson phase, depending 
on the dissipation parameter a = R q /Rs and the ratio 
Ej / Ec in the interesting parameter range. 

Our main conclusions can be summarized as follows: 
(i) Our detailed MC analysis unambiguously supports an 
earlier conjecture that in the zero temperature limit the 
superconductor-to-insulator phase transition always oc- 
curs at the dimensionless dissipation strength a = 1, in- 
dependently of the ratio Ej/Ec, (ii) at nonzero temper- 
atures this phase transition is substituted by a crossover, 
which position depends on the ratio Ej/Ec, as well as 
on temperature T and dissipation strength a, and (iii) 
at T — > this crossover line approaches the phase tran- 
sition line a = 1. These observations allow to fully rec- 
oncile the existing theoretical picture of the dissipative 
phase transition in a singje_resistively shunted JJ with 
the experimental results. E3E3 



II. QUANTUM DISSIPATIVE PHASE 
TRANSITION 



We proceed within the standard path integral formula- 
tion of the problem outlined in Rcf. ||[ The grand parti- 



tion function of the system "J J+shunt" can be expressed 
as a path integral over the Josephson phase <fi 



V</>exp(-S[</>]/h). 
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where S is the effective action 
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and P = h/hsT. The first and second terms in S[4>] 
account respectively for the charging and Josephson con- 
tributions. The third - nonlocal in time - term describes 
dissipation produced by an ohmic resistor. An additional 
dissipative contribution to the action, due to tunneling 
quasiparticles, is usually negligibly small in the interest- 
ing limit of energies/temperatures well below the super- 
conducting gap. Hence, in what follows this contribution 
will be ignored for the sake of simplicity. 

MC simulations with this effective action have been 
carried out by the standard discretization of the quantum 
paths into N imaginary-time slices.E3 In order to main- 
tain the accuracy of the calculation, as the temperature 
is lowered, the Trotter number N has been increased pro- 
portionally to the imaginary time (3. In our simulations 
we have taken N = 4f3Ec/h. This choice was proven to 
be sufficient in order to reach the necessary convergence 
of the calculated quantities: Repeating the calculations 
for some data points with larger N did not change our 
results. We have employed the Metropolis sampling to 
carry out PI MC simulations at temperatures down to 
ksT = Ec/125. A simulation run consists of successive 
MC steps, being updated all path-coordinates in each 
step. For each set of parameters (a, Ej/Ec, T), we 
generated 5 x 10 5 quantum paths for the calculation of 
ensemble-averaged values. More detaikr-ep this kind of 
MC simulations can be found elsewhere.E-3 

In order to quantify the quantum derealization of the 
phase we consider the mean-square fluctuations of the 
quantum paths 



((^) 2 ) = ((0(r)-0) 2 



(3) 



where <fi — J (j}(r)dT is the mean value of the phase 
for a given path. The quantity (|§|) was calculated by 
means of PI MC simulations with the aid of Eqs. (Til) and 
(!)■ 

In Fig. 1(a) we present ((S(f>) 2 } as a function of a 
for Ej /Ec = 2. Different symbols correspond to differ- 
ent temperatures; from top to bottom: ksT/Ec = 0.02, 
0.05, and 0.11. One observes that the phase derealiza- 
tion increases as the temperature is lowered. At a = 1 
one also observes a change of the slope in the data points. 
This change becomes more and more pronounced as the 
temperature is lowered, indicating the expected tendency 



2 



0.5 1 1.5 

a 



0.5 1 1.5 

a 



5 k 



to 

V 4 

3 
2 



-i 1 1 1 1 1 1 r - 



(b) 



Ej=0.75E c 



0.5 



1 

a 



1.5 



FIG. 1. Mean-square displacement ((Sep) 2 } as a func- 
tion of the dissipation strength a for (a) Ej = 2Ec 
and (b) Ej = 0.75Ec at three different temperatures: 
Squares, fcsT = O.llEc; circles, ksT = 0.05i?c; diamonds, 
k B T = 0.02E C . 



to a vertical line at a < 1 and T — > 0. Unfortunately the 
noise in the calculated values increases as T is reduced, 
especially for a < 1, where one expects an insulating 
regime (<j> delocalizcd). For a > 1, however, the quan- 
tum paths obtained in the MC simulations are merely 
confined to a single potential well, and the phase <f> is 
localized (superconducting regime). A similar picture is 
found for all other values of Ej/Eq at which our sim- 
ulations have been performed. The results obtained for 
Ej/Ec = 0.75 are presented in Fig. 1(b) for compari- 
son. In agreement with intuitive expectations, one finds 
that the phase fluctuations increase with decreasing ra- 
tio Ej/Ec- However, the phase transition point a = 1 
remains insensitive to this ratio in all cases. 

In Fig. 2 we show the average mean-square displace- 
ment ((S(f>) 2 ) versus a for several values of the ratio 
Ej I Eq at a temperature low as compared to the charg- 
ing energy: fc^T = 0.02Ec- Values of Ej/Ec increase 
from top to bottom: 0.2, 0.5, 1, and 2. As expected, 
for a given dissipation strength a, the phase becomes 
"more localized" as the ratio Ej/Ec increases, since the 
effective potential barrier for the phase increases as well. 



FIG. 2. Mean-square displacement ((8<t>) 2 ) as a function 
of the dissipation strength a for ksT = 0.02Ec- Results 
are given for several values of the ratio Ej/Ec- From top 
to bottom: Ej/Ec = 0.2, 0.5, 1, and 2. Error bars, if not 
shown, are of the order of the symbol size. 



Fig. 2 also demonstrates that stronger suppression of 
the phase fluctuations {(5(j)) 2 ) is always obtained when 
the dissipation strength a increases for a given value of 
Ej/Ec- Again in all cases we observe that the change of 
the slope in the data points for ((5<p) 2 ) occurs exactly at 
a = 1 , indicating the existence of the phase transition at 
this point. 

Further information can be obtained by investigating 
the dependence of ((6(f>) 2 ) on temperature. As expected, 
in our simulations we observe that for all values of E j / Ec 
and a, the quantum delocalization of <j> increases as the 
temperature is lowered. In Fig. 3(a) we have plotted 
the temperature dependence of ((S(f>) 2 ) for Ej/Eq = 2. 
Different symbols represent several a values, which in- 
crease from top (a = 0) to bottom (a = 1.2). At low 
T (Ec > lOfcsT), for all values of a one finds that 
((5<j)) 2 ) follows a power-law dependence on temperature, 
as shown by the dashed lines in Fig. 3(a). In particu- 
lar, in the dissipationless limit a = 0, it is reasonable to 
expect the quantum paths to be in a diffusive regime, so 
that ((5</>) 2 ) oc p. This is indeed confirmed by our calcu- 
lations. In the presence of dissipation the phase diffusion 
slows down, as indicated by a decrease in the slope of 
the dashed lines in Fig. 3(a). From this plot we find 
((S(f>) 2 } oc (P , with an exponent 7 smaller than unity 
for a > 0. Again, the localization phase transition at 
a = 1 is clearly observable in Fig. 3(a), since essentially 
no diffusion of the phase takes place at a > 1 and suffi- 
ciently large (3. Somcting similar occurs for other ratios 
Ej/Ec, as shown in Fig. 3(b) for Ej/E c = 0.75. The 
main difference with the previous case is that the low- 
temperature regime ((5</>) 2 ) oc (3 1 is reached at lower T 
{E c > 30k B T). 

In order to find the parameter 7, for each a 
we numerically evaluated the logarithmic derivative 
d\n((d(f)) 2 ) / d\n [3 = 7. The values 7 obtained in this 
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FIG. 3. Temperature dependence of ((S4>) 2 } for (a) 
Ej = 2Ec and (b) i5j = 0.75Ec at different values of a. 
From top to bottom: a = 0, 0.2, 0.4, 0.6, 0.8, 1, and 1.2. 
Dashed lines indicate the low-temperature trend of the data 
points. 



way are shown in Fig. 4 as a function of a for Ej /Ec = 
2. These results clearly indicate a linear dependence of 
7 on a of the form 7 = 1 — a. Thus, from our numerical 
analysis we can conclude that at sufficiently low temper- 
atures and a < 1 the phase diffusion is described by the 
formula 



(W>) 2 } oc 



(4) 



This dependence turns out to apply for all values Ej / Ec 
used in our simulations. From a linear fit to our data 
points in Fig. 4 we find the transition point at a — 
1.02 ±0.04. 

The prefactor in Eq. ([I|) depends on both Ej and 
Ec- In the limit Ej 3> Ec, onexan demonstrate with 
the aid of the instanton technique^ that this prefactor is 
proportional to A 1- ", where 
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FIG. 4. Low-temperature exponent 7 as a function of a 
for Ej — 2Ec- The dashed line is a linear fit to the data 
points. 
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is the bandwidth without dissipation. For Ej ~ Ec this 
equation does not apply anymore, but the qualitative 
trend remains the same: For a given value of Ec the 
prefactor increases monotonously with decreasing Ej\ cf. 
Figs. 3(a) and 3(b). 

To conclude this part of our analysis we emphasize 
again that - even though we present here our numerical 
results for the energy ratios Ej/Eq — 2 and 0.75 - sim- 
ilar behavior is observed for other values of the Joseph- 
son coupling energy. In particular, we have performed 
detailed MC simulations also for Ej/Ec = 0.25 and 3, 
and found essentially the same behavior as the one dis- 
cussed above. In all cases, from our numerical data we 
obtained unambiguous indications of the quantum local- 
ization phase transition at a = 1. 



III. EFFECT OF TEMPERATURE 

Now let us see how the above physical picture is mod- 
ified at nonzero T. Since any experiment is performed 
at a finite temperature, it is important to find out if the 
SIT can be observed under such conditions. 

To begin with, we recall the argumentta according to 
which superconductivity in Joscphson junctions can be 
observed even for a < 1, provided that the phase diffu- 
sion is slow enough, in the sense that the characteristic 
derealization time Tdei for the phase <f> exceeds the time 
of experiment. This situation can easily be achieved in 
the limit of very large Ej Ec- In that case the time 
Tdei is exponentially large, and no delocalization effects 
can be detected in a real experiment. This argument im- 
poses serious limitations on the observation of a SIT for 
sufficiently large values of Ej even at T — > 0. In Ref. 
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|l2| this argument was extended taking into account the. 
accuracy of the voltage measurements. The authorsEj 
argued that for their experiment l/rd e i should be com- 
pared with the quantity eV m i n rather than with the typ- 
ical experimental time, where Vmm is the minimum volt- 
age detectable in the experiment. They also noticed that 
further limitations can occur due to temperature effects. 

To explore the latter possibility we first notice that 
according to our result (^), at any nonzero tempera- 
ture quantum fluctuations do not spread the phase <f> 
to infinity even for a < 1. The phase does not simply 
have "enough time" to diffuse, and ((S(f>) 2 ) remains finite 
though increasing with the inverse temperature (3. Thus, 
at nonzero T and not very small Ej one might expect to 
observe a nonvanishing supercurrent even at small dissi- 
pation. This conclusion might appear paradoxical. One 
can argue that, if quantum fluctuations of the phase yield 
suppression of superconductivity already at T = (and 
a < f), at nonzero temperatures this suppression can 
only increase further, because of an additional effect of 
thermal fluctuations. 

In order to understand why this conclusion might not 
be quite correct, it is instructive to analyze the behavior 
of the (quasi-)charge variable!] canonically conjugate to 
the Josephson phase cj>. As discussed above, at T = 
and a < 1 the quasicharge is localized, i.e. Cooper pairs 
cannot tunnel across the junction due to Coulomb block- 
ade and, hence, the junction behaves as an insulator. At 
nonzero T this behavior persists as long as the temper- 
ature remains much smaller than the effective Coulomb 
gap for Cooper pairs. At Ej <C Ec this Coulomb gap 
is large (« Ec), while in the opposite limit Ej 3> Ec it 
turns out to be exponentially small:u 
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Here u p = y/8EjEc/h is the plasma frequency and A 
was defined in Eq. (ph. If the temperature becomes 
higher than the gap (™, fc^T > A r , Coulomb blockade 
for Cooper pairs (and, hence, the insulating behavior) 
is destroyed. This implies that the quasicharge Q gets 
strongly delocalized due to thermal effects. Because of 
the uncertainty relation 
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derealization of Q in turn restricts fluctuations of the 
canonically conjugate variable - the Josephson phase <f> 
- and as a result of that the superconducting behavior 
can be partially restored. The same scenario can be re- 
formulated in the phase space. One just needs to com- 
pare the typical inverse time during which the phase dif- 
fuses at a distance ~ 2tt with temperature. In the limit 
Ej 3> Ec this inverse time h/T& e \ ~ A r , and we ar- 
rive exactly at the same condition ksT ~ A r for the 
crossover line between "insulating" and "superconduct- 
ing" phases. Within logarithmic accuracy this crossover 
line agrees with the one found in Ref. 02. 
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FIG. 5. Dependence of the "Matsubara resistance" 
R(ui n ) oc |wn|(<^0)u„ on n — cj n p/2ir, for Ej/Ec = 2 and 
kBT — Ec/20. Different symbols correspond to different val- 
ues of the dissipation strength. From top to bottom: a — 0.6, 
0.7, 0.8, 0.9, and 1. 



Although the above consideration appears to be suffi- 
cient at high Josephson energies (Ej Ec), no quan- 
titative conclusion can yet be drawn for moderate cou- 
plings Ej ~ Ec- In order to study this parameter range 
one can apply a simple variational ansatzjj which leads 
to the self-consistency equations 

D = £jexp(-(<£ 2 ) tr /2) 
(0 2 )tr - fc s T^[CK/2e) 2 + aK|/27r + D}- 1 , (8) 



where u> n — 2im/(3 are the Matsubara frequencies. At 
T = these equations have a nonzero solution for D 
(which corresponds to superconductivity) only at a > I, 
whereas at nonzero temperature one can get a positive 
solution D > also for a < 1. By resolving these self- 
consistency equations at different T, one can qualitatively 
describe the above crossover for moderate values of Ej. 

This crossover can be studied quantitatively by PI MC 
simulations, using the effective action given above in Eq. 
(Q). With this purpose we have evaluated numerically 
the correlation function (0(r)^(O)) for different values 
of Ej/Ec, dissipation strength a, and temperature T. 
After Fourier transformation, this correlation function is 
directly related to the "Matsubara resistance" R(uj n ) = 
\ui n \ (4>4>) / Ae 2 , which yields the system resistance after 
analytic continuation to real frequencies. This numerical 
analytic continuation is a separate complicated problem, 
which will not be discussed here. Fortunately, however, 
this procedure is not needed in order to establish the 
position of the crossover line. 

Let us express the above correlation function in the 
form 
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are data points obtained from PI MC simulations. Dashed 
lines are guides to the eye. 



where Ej is the effective (renormalized) Josephson cou- 
pling energy. It is easy to see that in the low-frequency 
limit the behavior of the correlation function (^) is totally 
different, depending on whether Ej remains nonzero (su- 
perconductivity) or is fully suppressed by quantum fluc- 
tuations (insulating regime). In the first case, for suffi- 
ciently low frequencies the function (0) should inevitably 
decrease with u> n , while in the second case this function 
should increase and saturate at a finite value 2n/a in the 
limit tu n — > 0. Thus, by studying the behavior of the 
function (^]) it is possible to determine numerically from 
the PI MC simulations the position of the crossover line 
at different temperatures. 

In Fig. 5 we present our MC results for the func- 
tion |n|(#) w „ at Ej = 2E C , k B T = E c /20, and dif- 
ferent values of a. One observes that for sufficiently 
small n this function increases with decreasing n for 
a = 0.6 (upper curve), saturates for a = 0.7, and de- 
creases monotonously for a > 0.8 (three lower curves). 
Thus, by studying the small-n behavior of this function 
for different a values and Ej/Eq ratios, we arrive at the 
crossover line for the temperature UbT — Ec/20. Other 
temperatures are treated analogously. 

The resulting crossover lines are presented in Fig. 6 
for two temperatures: fc^T = Ec/5 and Ec/20. In 
full agreement with the above qualitative considerations, 
one observes that the position of the superconductor- 
insulator crossover line is shifted towards larger values 
Ej/Ec as the temperature is lowered. The same trend 
is observed for all other temperatures used in our simu- 
lations. Combining these results with those discussed in 
the previous section, one arrives at the conclusion that 
the crossover line should approach the phase transition 
line a = 1 in the limit T — ► 0. We would also like to 
point out that the position of the crossover line obtained 



within our MC analysis is fully consistent with that found 
experimentally in Refs. |l2||l3[ It appears, therefore, that 
deviations from the theoretical prediction for the phase 
boundary a = 1 observed in these exjieriments can be 
attributed to finite-temperature effects. c3 

In summary, the results of our MC simulations un- 
ambiguously demonstrate that the quantum (T = 0) 
superconductor-insulator phase transition in resistively 
shunted Josephson junctions occurs at the value of the 
shunt resistance R$ = R q , irrespective of the ratio 
Ej/Ec- For a < 1, quantum diffusion of the Joseph- 
son phase 4> yields a simple scaling dependence {(54>) 2 ) oc 
Finite-temperature effects turn the phase transi- 
tion into a crossover, which position depends on the ratio 
Ej/Ec, as well as on the dissipation strength a, and on 
temperature. Our resulls_are fully consistent with recent 
experimental findings £3113 
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